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Abstract. We derived the properties of V605 Aql before the final helium flash pulse by studying its surrounding PN A58. 
Photoionizing models of our spectral data together with a new distance estimate and a closer look at the recombination 
timescales lead to a consistent model. Comparing our findings with the only hydrogen-poor twin, namely Sakurai's Object, 
we conclude that these born-again objects have normal PNe core masses. We are able to prove V605 Aql indeed to be, similar 
to V4332 Sgr, a very late thermal pulse object and to put constraints for the evolutionary time scales for the transition back to 
the AGB. 
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1. Introduction 

Calculations done in the early 1980s demonstrated that lumi- 
nous central stars of PNe could be undergoing a final helium 
shell flash before the final cooling of the white dwarf. It was 
reviewed by (Iben & MacDonald 1995 1 that this last thermal 
pulse brightens a white dwarf to AGB luminosity and results 
in a sudden change in abundances to heavier elements as the 
residual hydrogen envelope is consumed in the helium-burning 
convective shell. Back at the AGB the central star of the PN 
(CSPN) pursues again its evolutionary path to become a white 
dwarf, but this time the main energy source is helium burn- 
ing. This evolutionary behavior is called the born-again sce- 
nario. Born-again PNe helium-flash objects are characterized 
by highly processed hydrogen-deficient nebular material and 
dust and are of fundamental importance for the understand- 
ing of the late stages of stellar evolution. Due to the short 
timescales of the transition only few born again objects were 
discovered therein. At the moment there are actually only three 
objects - Sakurai's Object (V4334 Sgr), V605 Aql and FG Sge 
-, where we can observe more than a final H-deficient PN. This 
makes individual detailed studies of these objects necessary, 
especially as they are also discussed as possible progenitors of 
[WC]-CSPNe, PG1195 type CSPNe and RCrB stars. For a re- 
view see Clayton ( 1200 It and references therein. 

Sakurai's Object (V4334 Sgr) is the proto-typical final he- 
lium shell flash object and had its flash in 1996 (Diirbeck & 



Bene tti [1996). The old PN of V4334 Sgr has been found and 
its spectrum reveals it to be an ordinary PN (Pollacco" 1999V 

V605 Aql is often refeiTed to as the 'older twin of Sakurai's 
object', as many similarities exist. It was discovered in 1919 as 
'Nova Aquilae No.4' and later assigned V605 Aql. Only in the 
1980s was it recognized to be a final helium flash object and 
thus it was observed less than Sakurai's object or FG Sge. It 
had its flash already in 1917, and therefore serves as pathfinder 
for future events in V4334 Sgr, which is still frozen in its reborn 
AGB resp. post-AGB state. A detailed historical review of the 
evolution of V605 Aql can be found in Clayton & de Marco 
Cl997j and the photometric similarities between V605 Aql and 
V4334 Sgr are discussed in Diirbeck et al. 
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Different theoretical models exist that should describe the 
born-again phenomenon, but there is not yet a consent about 
the crucial parameters like convective mixing, opacities at high 
ionization energies, the importance of convective overshoot 
or surface cooling by means of stellar wind (e.g. Lawlor & 
MacDonald 2003 Herwig2001l. 

We present here a detailed study of V605 Aql and compare 
our results to those found for V4334 Sgr 



2. Observations and Data Reduction 

We obtained spectra using the ESO NTT telescope at LaSilla 
(August 4* to 6* 2002) with the multi mode instrument EMMI 
mounted. Also different narrow-band filter images were ob- 
tained there. The calibration (bias, flatfield, wavelength cali- 
bration and response curve for the spectra; bias and skyflats for 
the images) was done using usual procedures in MIDAS. The 
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resolutions of the individual grisms used for the spectra can be 
found in Tabled 

Table 1. Grisms used for spectroscopy 



Grism 



Resolution 



grism #2 
grism #3 
grism #6 



0.35 nm/pixel 
0.29 nm/pixel 
0.15 nm/pixel 



390 to 980 nm 
390 to 850 nm 
500 to 880 nm 









Fig. 1. These [OIII] (left) and [Nil] (right) narrow-band images 
of PN A58 are not rotated to prevent information loss. The po- 
sition angle to a north-east orientation is 132.381. The edge en- 
hancement in [Nil] gives the impression of an interaction with 
the ISM. 

Narrow-band direct images (see Fig.^ show a strong de- 
pendence of the shape of PN ASS on the excitation level. The 
[OIII] image of the PN looks rather roundish, whereas in H/S 
and in [Nil] a western rim is clearly visible. In all three im- 
ages the nebula seems to consist of three parts; an inner central 
part and two 'polar caps', which is most obvious in the [OIII] 
image. According to the manual the EMMI CCD chip has a 
resolution of 0. "333 per pixel, but the seeing was not sufficient 
to derive the size of the V605 Aql knot as done by Hinkle et 
al. (J2001). On the basis of the CCD images our estimate for the 
size of the old PN A58 is 50" for the longitudinal axis and 38" 
for the secondary axis. 

The edge enhancement suggests an interaction with the in- 
terstellar medium (ISM). Hence we determined the physical 
properties of the old nebula Abell 58 out of the inner part, 
nearer to V605 Aql. 

The interstellar extinction had been derived from the 
Balmer decrement. Different regions show an E(B-V) with ex- 
trema of 0'."35 and 0'."55. A weighted average of EiB-V) = 0'?5 
was adopted for the models. 

3. Distance Determination 

There were several distance estimates for V605 Aql before, 
giving values between 2.7 and 6.0 kpc as reviewed by Clayton 
& de Marco ( 1997 1. who recommended a distance of 3.5 kpc. 
We tried to determine the distance by the use of radio ob- 
servations. The VLA Sky Survey by (Condon & Kaplan 1998^ 
gives a S - 2.7 mJy flux-density peak for PN A58 at 1.4 GHz 



(according to a wavelength of 21 cm), with an 0.8 mJy RMS 
uncertainty. Usually distances are calculated by the 6 cm radio 
fluxes (according to a frequency of 4885 MHz), but assuming 
free-free radiation with 



/ = /()V 

we translated this value into 



-0.1 



* 6cm 



0.0024 + 0.0007 Jy . 



Now the method by Cahn et al. J 19921 1 can be applied. The 
optical thickness parameter t is given by 



T = log 



402 



* 6cm 



with S(,cm in Jansky and 6 in arcseconds. A mean value for the 
angular radius of the PN of 22" resulted in t = 5.91+0.12. For 
T > 3.13 the ionized mass yU of the PN can be described with 

log(/i) = -0.87 . 

Using these equations, the radius R and the distance D of the 
PN can be determined from the ionized mass by 

log(/.).l||> -1.306+1 =-1.654+1 



and 



D = 206.265 - . 



So the solutions using (Condon & Kaplan [T998> and (Cahn et 
al. [T992l are 

R = 0.34 + 0.02 pc 



and 



D = 3.15+j;:?^ kpc. 



Alternatively the method of van de Steene & Zijlstra J1995> 
uses the radio continuum brightness temperature 



Tb 



^ ^ 6cm 1 o Af\f\ ^cm ^ ^^ , , -, f\-2 



2nkv^ 6*2 



18400- 



02 



9.05 X 10" 



From this temperature, which is a measure for the brightness of 
the nebula at radio frequencies, together with their fitted linear 
regression line through a sample of bulge PN they specify a 
radius R: 

log(R) = -0.35(+0.02) logiTh) - 0.51 (+0.05) 

leading to R - 0.71 pc and D - 6.6 kpc. This approximation is 
generally a factor of 2 too big. 

Normally the best and most faithful results (for a discussion 
see Schmeja & Kimeswenger 2002 1 are derived by the method 
of Schneider & Buckley ( 1996. 1. Their equation is: 

log(D) = - log(0) - 0.0261 log2(/) - 0.299 log(/) + 1.198 , 

where the surface brightness / = S(,cmlifi^'^) is measured in 
mJy arcsec"2 xhis leads to D = 3.08+y[2 kpc. 

For all further calculations we therefore assumed a distance 
of D = 3.1 kpc. This means that the radius R of the PN A58 
equals 0.33 pc. 
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4. Photoionization Models 

On the basis of the line fluxes presented in Table|2lthe old PN 
A58 was modelled with CLOUDY (Ferland 1996 1. The final 
goal was not to reproduce the ionization and thermal equilibria 
within the nebula in detail but to qualify the stellar continuum 
before the helium flash pulse. The low density PN sees still its 
old central star at the time, when it was sufficiently hot to ionize 
it. Thus a large grid of photoionizing models was calculated, 
giving us the position of V605 Aql in the Hertzsprung-Russell 
diagram (HRD) during its hidden past. 

Table 2. Dereddened Une fluxes for the PN A58 

Line ID A [A] Intensity [U/S = 100] 

27 
100 
286 
804 
4.5 ±1.0 

18 

92 
286 
280 ±20.0 

42 

31 

13 



Usually blackbody radiation is used for such calculations 
of other objects (e.g. van de Steene '1995' van Hoof & van 
de Steene [T9991 Kerber et al. 1999 Niirnberger et al. 2001 
Cuisinier et al. l2002> . Armsdorfer et al. J2002> and Armsdorfer 
& Kimeswenger (> 2001> showed that the use of real stellar at- 
mospheres changes especially the helium lines. Here we thus 
used NLTE central star models of Rauch ( flWTI |2003 1. The 
same grid of model atmospheres was used also by Pollacco 
(|TH55) for the twin object V4334 Sgr. 

The temperature of the central star of the PN, its luminos- 
ity, the gas density and the filling factor of the gas were varied. 
While the range for the filling factor was selected to fit to the 
typical findings for evolved nebulae (Armsdorfer et al. 2002 1, 
the gas density range of 100 to 300 e^cm"^ was taken from 
the ratio of the [SII] doublet as in Osterbrock (1974 i. The dis- 
tance of 3.1 kpc was adopted from section |51 leading with the 
average angular radius of 22" to the linear radius of the shell. 

The strength of the helium lines is nearly independent from 
the CSPN luminosity. Thus this fixes the abundance of this el- 
ement clearly. The variations of the parameters were chosen 
to get a optimized intersection of the solutions for the Helssve, 
the [OIII]4948+5007 and the [NII]6548+6583 lines at the level of 
^bserved/^iodei = LO. These Solutions are not completely in- 
dependent from the chosen abundance. The [Nil] 5755 line then 
should intersect, if the resulting electron temperature is suited, 
there as well. This fixes the CSPN temperature and therefore 
the whole solution. The best fit model is shown in Fig.|2] 

The placement of the derived CSPN luminosity and tem- 
perature together with the theory of stellar post-AGB evolution 
by Blocker (1995) yields the mass of the final white dwarf and 
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Fig. 2. The best fit model shows a good intersection of the so- 
lutions for the Helium, Oxygen and Nitrogen lines at the level 
^bserved/4iodei - 1-0. The resulting model is given in Tabled 
Table 3. Parameters of the best fit model. 



Parameter 



Value 



distance 

CSPN: 
luminosity 
effective temperature 
age from post-AGB 
white dwarf mass 

PN gas shell: 
average radius 

thickness 
filling factor 
electron density 
abundances: 

helium 

nitrogen 

oxygen 



3.1 kpc (see section|3j 

L = 325±100Lo 
Teff = 121.5 ±6kK 
t = 8300!^^°" years 
M = 0.605 ± 0.02 Me 

angular 22" 

lineal- ?-= 1.02 10'* m 

A,/r = 0.2 

6 = 0.05 (adopted) 

n, = 200 ± 100 e-cm-3 

[He] = -0.83 
[N] = -4.00 
[O] = -3.40 



the transition time t. The position of the CSPN in the HRD 
is shown in Fig.|3] The resulting model (see Table |3} is com- 
pletely self consistent. In addition it is also supported by the 
velocity measurements by Pollacco et al. J1992> . who derived 
for the expansion velocity of A58 v^^p = 31+4 km s~' , giving 
with the dynamical age of 8200 years a dynamical radius of PN 
A58of 

R = 0.26 ±JJ:Jg pc , 

which matches extremely well our radius of R = 0.33 pc as 
derived in section |3 This is an especially good conjunction, if 
we take into account that the dynamical solution is usually a 
factor of 2 or 3 off target. 

5. Recombination Timescales in A58 

There is definitely Hell4685 in the spectrum of A58. This 
fact provides together with the knowledge of recombination 
timescales an estimate for the time between the ending of the 
UV radiation of the central star and the nova burst (position H 
or E inFig.|3}- 
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Fig. 3. The position of the CSPN in the HRD at the moment 
of the helium flash. The evolutionary tracks for the post-AGB 
evolution for a 0.605, 0.625 and 0.696 M© white dwarf from 
Blocker ( 1995 1 are plotted. For comparison reasons also the 
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but enlarges the [OIII] recombination time-scales. For a more 



results of PoUacco S199% and of Kerber (IT999b for Sakurai's 
Object are shown. Pollacco's error bars show the observational 
incertainties, he declares log{L/LQ) < 3 to be consistent with 
the white dwarf cooling tracks. Kerber assumed a relative error 
of 10 percent in the intensities, but admits that the overall error 
bars are certainly larger. The grey line underneath represents 
the evolutionary model of Lawlor & MacDonald J2003> . 

Using the same dielectric recombination coefficients a as 
PoUacco <1999> we calculated a table for the recombination 
time-scales for electron densities n^ of 100 and 200 e^cm"^ 
(see Table|4li. 

Table 4. Recombination time-scale for atomic states at T^ = 
10.000 K for electi-on densities rig of 100 and 200 e cm ^ An 
increase of the electron density to n^ - 300 e cm"^ only lowers 
the upper limits. The terms 'r' and 'd' refer to the radiative and 
dielectronic components in the total recombination coefficient 
't'. 



State T, 


(n, = 100) [yr] r 


(n, = 200) [>t] 


H" 


760 


380 


He" 


1500 


750 


He+ 


206 


103 


02+(r) 


31 


16 


02+(d) 


28 


14 


02+(t) 


19 


9 


O^r) 


155 


78 


O^d) 


190 


95 


0+(t) 


85 


43 


NMr) 


139 


70 


N+(d) 


155 


78 


N+(t) 


73 


37 



It can be seen that the recombination times for [OIII] are 
relatively low, but this table is only valid for single atomic sys- 
tems. We have to take into account the Bowen Hell Lyman a - 
[OIII] resonance-fluorescence mechanism within the PN. This 
mechanism slightly decreases the Hell recombination times 



detailed study refer to Osterbrock (I1974> . The emphasis lies 
upon Hell as most common element, but also Nil is important 
as it is not influenced by the Bowen Hell - [OIII] resonance- 
fluorescence mechanism. 

Since we have to take the minimum of all the derived time- 
scales for a specific electron density as an upper limit for the 
lifetime of a nebula in this state, we can use a weighted average 
of 100 + 50 years as a good estimate for the upper limit. This 
also justifies the use of photoionizing models with CLOUDY. 

6. Discussion and Conclusion 

The parameters used to calculate stellar evolution become 
crucial when we consider such subtleties as born-again objects. 
Hence the physical properties of born-agains in all stages pro- 
vide a hard constraint on valid evolution theories. In general 
two different scenarios are discussed for the born again event: 

- The late thermal pulse (LTP) where the helium flash occurs 
while the star is still H-shell burning. The CSPN has stifl 
a luminosity above a few thousands L©. H-ingestion to the 
flash convection zone causes - according to current models 
- a slow evolution back to the AGB lasting at least several 



hundreds of years (Blocker & Schonberner l200 1 l l . 
- The very late thermal pulse (VLTP) where the CSPN was 
akeady on the white dwarf cooling track, i.e., after the ces- 
sation of H burning. In this scenario the stellar luminos- 
ity decreases slowly from the flash position (position A in 
Fig. 13 to a very low luminosity (called B in Lawlor & 
MacDonald 2003 1. The timescale for this transition is dis- 
cussed very controversially. While Lawlor & MacDonald 
have transition timescales of up to several thousands of 
years, Iben et al. (1983"l and Iben & MacDonald (1995"l give 
a few tens of years, Herwig et al. (1999) and Herwig ( 2 0011 
I2003> give a timescale of less than one year for a 0.604 Mq 
CSPN. Older models give such a fast evolution only for 
high mass CSPNe. 



For a more detailed discussion see Blocker J2001t and refer- 
ences therein. 

Based on recombination timescales we conclude that the 
time between the end of the UV radiation of the central star 
(tfiash - the time of the late helium flash) and the observed 
bright Nova burst of V605 Aql in the 1920's (= f„„„,) has to 
be very short. Modelling the suiTounding old PN A58 with 
CLOUDY, we obtain that the central star had a mass near 
0.6 Mq (using the tracks of Blocker (1995) as well as those 
of Herwig et al. (199 9 ) and Herwig ( 2001 ) ) and a temperature 
of 120000 Kelvin before the final helium flash. The model of 
Herwig ( 2001") with lower masses is neither consistent with the 
photometric data of V4332 Sgr, as pointed out there already, 
nor with our model of the pre-outburst CSPN of V605 Aql. 
According to Blocker (1995") the time since leaving the post- 
AGB resulted in 8200 years. The distance of 3. 1 + 0.2 kpc from 
radio observations goes well with the constraints by the pho- 
toionization model and also matches the angular size for the 
derived age. Our observations obtained in 2002 (- tabs) and 
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thus about 85 years after t„oya and the finding that the age of 
the PN (fobs - tfiash) has to be around 100 years, lead to the 
conclusion that the transition time to either point H or E has 
to be below 20 years - taking into account all types of uncer- 
tainties. This clearly excludes a LTP and accordingly is consis- 
tent with the CSPN position in the HRD. Our findings clearly 
- both in luminosity and especially in temperature - exclude 
the possibility that the PN was excited in position D or F (see 
Fig. O. In addition, test calculations with CLOUDY also ex- 
clude the possibility that the PN was excited additionally in 
between by a 40 000 < T(-|5pj>j < 70 000 K source with a 
luminosity more than IOOOLq. This would have imprinted by 
special emission line ratios of [Nil] and [Sll]. Also, although 
checked by Fuhrmann (198 1 ) on more than 430 Sonneberg and 
Heidelberg patrol plates taken after 1920 and by Diirbeck et al. 
<2002i) on the Harvard patrol plates, no brightening has been 
detected since the 1917 outburst. This implies that the source 
is either in its first return to the AGB (position E), or not visi- 
ble due to a thick dust shell during the transition E-F-H, or this 
second loop does not really occur. Whilst Iben et al. J1983> do 
not show such a loop at all, the transition time for the model 
of Lawlor & MacDonald (2003) gives timescales of only a 
few tens of years for the heating process after the first "Nova 
event" for the transition E-F. The timescale in Herwig (2003| 
is even shorter. Then it stays near F for a long time at high lu- 
minosity and high temperature. This environment would have 
destroyed a dense dust shell built during the short evolution 
from leaving E towards F (Koller^OOQ, KoUer & Kimeswenger 
I2001all2001bl Kimeswenger & Koller l2002l l and consequently 
V605 Aql should have been observable as a blue star on the 
POSS plates. This wasn't found neither on the first genera- 
tion POSS plates in the 1950s nor on the second generation 
in the 1980s and 1990s. We thus conclude that the object is 
at its first return to the AGB and the timescale for the transi- 
tion has to be as short as given by Herwig ( 1200 1> and not as 
long as given by Lawlor & MacDonald (2003). Furthermore 
V605 Aql should have stayed at least several decades at the 
AGB before reheating to the temperature seen nowadays in 



the CIV lines (Kimeswenger 120031 Lechner & Kimeswenger 
12003 1. This gives somewhat longer timescales than those pre- 
dicted by Lawlor & MacDonald UM^ or by Herwig ( I2003> . 

The mass of the progenitor CSPN also was discussed 
controversially. Diirbeck et al. (.1997) derive >0.7 M© using 
Paczynski's core-mass relationship. Diirbeck et al. ( I2000> ar- 
gue, using the high luminosity and the evolution time span, 
that a mass of 0.8 Mq should be a lower limit. However the 
outcomes here for V605 Aql and those of PoUacco 119991 1 and 
Kerber (T^W) (see Fig. place V605 Aql and V4334 Sgr 
clearly in the range of normal PNe core masses at or just be- 
low 0.6 Mq. Kerber ( I1999t calculated his models with the pho- 
toionization code GWYN and used only blackbody radiation. 
Pollacco (1999) used CLOUDY with different but compara- 
ble stellar atmospheres. It can be assumed that their spectra 
are from different regions in Sakurai's Object. In Fig.|3we see 
also beautifully the dependence on the distance with the aid of 
Kerber's results for 1.5 kpc and 5.5 kpc for Sakurai's Object. 
The mass of 0.6 M© also corresponds extremely well with the 
dynamic age of the old PN. A CSPN with a mass above 0.7- 



0.8 M0 would evolve much faster than the observed PN. So we 
conclude that the CSPN mass is very similar to that of Sakurai's 
object found by Herwig ( 1200 1> . 
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